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Abstract 

The concept of intermediate tail dependence is useful if one wants 
to quantify the degree of positive dependence in the tails when there is 
no strong evidence of presence of the usual tail dependence. We first 
review existing studies on intermediate tail dependence, and then we 
report new results to supplement the review. Intermediate tail de- 
pendence for elliptical, extreme value and Archimedean copulas are 
reviewed and further studied, respectively. For Archimedean copulas, 
we not only consider the frailty model but also the recently stud- 
ied scale mixture model; for the latter, conditions leading to upper 
intermediate tail dependence are presented, and it provides a useful 
way to simulate copulas with desirable intermediate tail dependence 
structures. 

Key words and phrases: Tail order, elliptical copula, extreme value 
copula, Archimedean copula, frailty model, resilience model, scale 
mixture model. 

1 Introduction 

For applications in many areas such as environmetrics, actuarial science and 
quantitative finance, in addition to prudent examinations of univariate mar- 
gins, careful modeling of various dependence patterns in corresponding distri- 
butional tails is often very important. For statistical modeling of dependence 
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structures between random variables, a very useful approach is to employ a 
copula function to combine the univariate margins together to get the joint 
distribution. A copula C : [0, l] d — > [0, 1] for a rf-dimensional random vec- 



tor can be defined as C(ui : 



.u d 



F(F^( Ul ) 



p-h 



(ud)), where F 

is the joint cumulative distribution function (cdf), Fj is the univariate cdf 
for the zth margin, and F~ x is the generalized inverse function defined as 
Ff (u) = inf{x : Fi(x) > u}. To avoid technical complexity, throughout the 
article, the univariate cdfs Fj's are assumed to be supported on [0, oo), and 
Fj's are continuous and thus the copula C is unique due to Sklar's theorem 
(Sklar, 1959). We refer the readers to Joe (1997) and Nelsen (2006) for refer- 
ences of copulas. Moreover, all distribution functions and density functions 
are assumed to be ultimately monotone to the left and right endpoints; this 
condition is very mild and satisfied by all the commonly used distributions. 

The so-called tail dependence parameters (also called tail dependence 
coefficients or tail dependence index) have been studied as a summary quan- 
tity to capture the degree of tail dependence. Let A be the upper tail de- 
pendence parameter, then A := lim u _ 5>0 + C{u, . . . , u)/u (provided that the 
limit exists), where C is the survival copula of C\ that is, C(u\, . . . ,Ud) '■= 
C(l — Ui, . . . , 1 — Ud), and C is the survival function of C and defined as 
C(«i, . . . ,u d ) := 1 + X^/c{i,...,4( _1 ) |/|c V( u *' i e -0> where Cj is the copula 
for the /-margin. 

The concept of intermediate tail dependence arises when one wants to 
quantify the strength of such dependence in the tails but the usual tail de- 
pendence parameter A = 0. In this case, one needs to find another quantity 
to capture the strength of dependence in the tails. 

In Hua and Joe (2011b), the concept of tail order (k) is suggested as a 
quantity to capture the leading information of dependence in the tails when 
A = 0. When we look at the decay of a copula function along the diagonal 
as a function in it, then a mild assumption is that C(u, . . . ,u) ~ u K £(u) as 
u — > + , where i is a slowly varying function (Bingham et al. , 1987) and 
the notation g ~ h means that the functions g and h are asymptotically 
equivalent; that is, lim t _j. to g(t)/h(t) = 1 with t being the corresponding 
limiting point that is usually or oo. The leading parameter k is referred to 
as the lower tail order of copula C. In parallel, the upper tail order of copula 
C can be defined as the n that satisfies that C(u, . . . ,u) ~ u K £(u),u — > + . 
Clearly, k = 1 corresponds to the usual tail dependence and then the limit of 
the slowly varying function £(u) is used as the quantity to capture the degree 
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of tail dependence. 

If copulas C\ and C 2 have lower tail orders K\ and k 2 respectively and 
Ki < K 2 , then C\(u) > C 2 (u) for all u in a neighborhood of 0. Similarly, 
this holds for the upper tail orders in terms of survival functions. That is, 
the tail order for comparing copulas implies a form of multivariate stochastic 
order in the joint lower or upper tails. 



To the best of our knowledge, Ledford and Tawn (1996) is the first paper 



that employs a regularly varying function to study the weaker dependence 
in the tails with the tail dependence parameter A = 0. More specifically, for 
a bivariate random vector (Xi,X 2 ) T , where X\ and X 2 are unit Frechet dis- 
tributed with cdf Fi(x) = e~ x l x x > 0, i = 1, 2, and are nonnegatively associ- 
ated, assume P[Xi > r, X 2 > r] ~ £(r)r -1 /' 7 ,r — > 00, where 1/2 < 77 < 1. It 
can be verified that the tail order k corresponds to 1/rj of Ledford and Tawn's 



representation, and r\ is called the residual dependence index in Hashorva 



(2010) and references therein. A lot of research has been done following this 



Tawn 


(1996 


1997 


); 



direction. We refer to Coles et al. (1999); Heffernan (2000); Ledford and 



this idea. 

Although the concept of tail order is defined with respect to a copula, 
which is a more intuitive way, the k itself or some functional forms of k 
simply describe the relative speed of decay of the joint tail probability to 
certain functional forms of the tail probability of one of the standardized 
margins. How to standardize the margins and how to choose the functional 
forms depend on how to make such relative speed of decay meaningful and 
be able to capture the leading information of dependence in the tails. 

Another notion that is close to the concept of tail order is a tail depen- 
dence measure in the sense of Coles et al. (1999), in which an upper tail 



dependence measure for a bivariate copula is defined as 

X ■= lim 2 log(w)/log(C(l — u, 1 — u)) — 1. 

So, x — 2/ K u ~ 1- Note that if C(u, . . . ,u) ~ u K £(u) as u — )■ + , letting 
C'(u, ...,u) := d(C(u, . . . ,u))/du and [u K £(u)] f := d(u K £(u))/du, then by the 



l'Hopital's rule and the Monotone density theorem (Bingham et al. , 1987) 
\og(C(u,...,u)) 



lim 

u->0+ 



log(w) 



, C'(u,...,u) u\uH(u)Y 

hm — r — X 



u ^o+ [u K £(u)}' 
ku r £(u) 



C(u,...,u) 



lim ru 

u^0+ C (u, . . . , u) 



K. 



(1) 



3 



Therefore, in some cases it may be easier to obtain the lower tail order by 
applying Q; for example, it will be used in the proof of Proposition [3] for a 
bivariate elliptical copula. Similarly, the upper tail order can be calculated 
as K V = lim u ^ +[log(C(l 1 - u))]/pog(it)]. 

In what follows, we will introduce the concepts of tail order and inter- 
mediate tail dependence in Section [2j Some detailed results for extreme 
value, elliptical and Archimedean copulas are presented in Sections |3j [4] and 
|5j respectively. In particular, in Section [5j we will study intermediate tail 
dependence through two different stochastic representations of Archimedean 
copulas: the frailty model and the scale mixture model. The study of upper 
intermediate tail dependence of Archimedean copulas derived from the scale 
mixture model is new in the literature. References are given for existing re- 
sults, and proofs are only provided for new results. Section [6] will conclude 
the article. 



2 Tail order and intermediate tail dependence 



The theory of regular variation will be applied throughout the article. We 



refer the reader to Bingham et al. (1987), Resnick (1987, 2007), de Haan and 



Ferreira (2006) and Geluk and de Haan (1987) for references. A measurable 



function g : IR + — > IR + is regularly varying at oo with index a^O (written g G 
RV a ) if for any t > 0, \im. x ^. 00 (g(xt) / g(x)) = t a . If the above equation holds 
with a = for any t > 0, then g is said to be slowly varying at oo and written 
as g G RVV For regularly varying at 0, that is, \im x ^ +(g(xt) / g(x)) = t a 
for any t > 0, the notation is g G RVq,(0 + ), and slow variation of £ at is 
written as £ G RV (0 + ). For any g G RV a ,a G E, there exists an £ G RV 
such that g(x) = x a £(x). 

Other notation: a bold letter is used to represent a transposed vector, 
e.g., x := (xi, . . . , Xd), and (uld) represents (u, . . . , u) with d components of 
it's and (uw) represents (uwi, . . . ,uwd). We use (wj) := (wi,i G I) when / 
is a subset of {1, . . . , d}, and sometimes (u>i, . . . , wo) is abbreviated as (w d ). 

Definition 1 Suppose C is a (^-dimensional copula. If there exists a real 
constant Kl(C) > and £ G RVq(0 + ) such that 



C(ul c 



u 



u 



u — y 



then we refer to kl{C) as the lower tail order of C and refer to A^(C) = 
lim M _ s>0 + C-{ u ) as the lower tail order parameter, provided that the limit exists. 
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Similarly, the upper tail order is defined as kjj(C) that satisfies C((l— ii)l<f) ~ 
u Ku ( c H(u),u + , with the upper tail order parameter being A[/(C) = 
lim u _ s>0 + ^(w), provided that the limit exists. 

The notion of tail order is especially useful for bivariate copulas or for 
multivariate permutation symmetric (i.e., exchangeable) copulas. Otherwise, 
the study of tail order involves more technical issues. All the copulas studied 
in this article are assumed to be permutation symmetric, in order to illus- 
trate the main ideas and key results without involving too much technical 
discussion. Under such assumptions, intermediate tail dependence simply 
means that the corresponding tail order k satisfies 1 < k < d. 

Definition 2 Suppose C is a d-dimensional copula and C(uld) ~ u K £(u), 
u — > + for some £ G RV (0 + ). The lower tail order function b : R+ — y R + is 
defined as 

K«iO,»)= ^ 

u^o+ u K £(u) 

provided that the limit function exists. In parallel, if C((l — u)ld) ~ u K £(u), 
u — y 0^ for some £ G RVo(0 + ), the upper tail order function b* : R^ — y R + is 
defined as _ 

m-s>o+ u K £(u) 
provided that the limit function exists. If £{u) — y h =fi 0, then hb(w; C, 1) 



and hb*(w; C, 1) become the tail dependence functions proposed in Joe et al. 



( 2010[ ). That is, the definition of the tail order function here absorbs a 



constant into £ so that 6(1^; C, k) = 1. 

The following are some elementary properties of the lower and upper tail 
order functions b and b* . Obvious properties of tail order for C are the 
following: k l (C) = Ku{C), Ku(C) = k, l (C), b(w;C,K) = b*(w;C,K) and 
b*(w; C, k) = b(w; C, k). 

Proposition 1 A lower tail order function b(w) = b(w; C, k) has the fol- 
lowing properties: (1) b(w) = if there exists an i e {1, . . . , d} with Wi = 0; 
(2) b{w) is increasing in Wi,i G {l,...,d}; (3) b(w) is positive homoge- 
neous of order k; that is, for any fixed t > 0, b(tw) = t K b(w); (4) if 
b(w) is partially differentiable with respect to each Wi on (0, +oo) ; then 
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Other properties of tail order and tail order functions can be found in 
Hua and Joe (2011b). Due to the limitation of space, in what follows, we 



will focus on the study of intermediate tail dependence for three important 
copula families. 



3 Extreme value copula 



If a copula C satisfies C(u\, 



C*(mi, ...,u d ) for any (u u ...,u d ) e 



[0, 1] and t > 0, then we refer to C as an extreme value copula. For any 
extreme value copula C, there exists a function A : [0, oo) d — > [0, oo) such 
that 



C(ui, ...,u d )= exp{-A(-logwi, . . . ,-\ogu d )}, 
where A is convex, positive homogeneous of order 1 and satisfies 
max(xi, ... ,x d ) < A(xi, . . . , x d ) < xi H h x d . 



(2) 



We refer to |Pickands| ( |1981[ ), |Tawn| ([1990]) and Chapter 6 of |Joej (|1997j) for 

references. 



For the lowe r tail of an extreme value copula C, it can be verified (Hua 

u^ 1 '"-' 1 '. That is, for any extreme value 



and Joe, 2011b) that C(u, . . . ,u 
copula C, the lower tail order is kl{C) = A(l,...,l), and there is inter- 
mediate lower tail dependence except for the boundary cases, such as the 
independence copula and the comonotonicity copula, where A(l, . . . , 1) = d 
and 1, respectively. 

For the upper tail of an extreme value copula C, write the survival copula 
as C(ul d ) := C((l - u)l d ) = 1 + E ^c { i,..., rf} (-l) |7| ^((l ~ Since 



each /-margin of C with 2 < |/| < d is also an extreme value copula (Tawn 



1990), let Ai be the corresponding function in the sense of (|2| for the extreme 
value copula C7, then, 



C(u,...,u) 



C{l-u,..., 
l-d + du + 
1 — d + du + 



1 -u) 



1 + 



E 

0^IC{1,.. 



J2 (-l) w C z (l-u m ) 

IC{l,...,d},\I\>2 

E ( 

/C{l,...,d},|7|>2 



u 
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By the l'Hopital's rule, as u — > + , 



c{u,...,v)~u\d- J2 (- 1 ) m ^/( 1 m) = :uX - 

\ IC{l,...,d},\I\>2 J 

So if A 7^ 0, then C has usual upper tail dependence with tail dependence 
parameter \ = d — Y2icn 4 |/|>2( — If A = 0, then it is unclear 
if an extreme value copula can have upper intermediate tail dependence or 
not; certain structures of those Aj functions are needed in this regard. 



4 Elliptical copula 

Since a copula is invariant to a strict increasing transformation on margins, 
for the study of elliptical copula, we may omit the location and scale param- 
eters of joint elliptical distributions. Intermediate tail dependence depends 
just on the radial random variable, and the condition on the radial random 
variable can be seen from the bivariate case, so the main result in this section 
is bivariate, and it can be easily extended to exchangeable multivariate ellip- 
tical copulas. Now consider the following representation: let X := (Xi,X 2 ) 
be an elliptical random vector such that 

X = RAU, (3) 

where the radial random variable R > is independent of U, U is an bi- 
variate random vector uniformly distributed on the surface of the unit hy- 
persphere {z e IR 2 |,z T z = 1}, A is a 2 x 2 matrix such that AA T = S where 
the entries of S are En = S 2 2 = 1 and S 12 = E 2 i = g with — 1 < g < 1, e.g., 

A = ( r— — k ] . For such an elliptical distribution, the margins have the 
yg^/l-gy 

same cdf assumed to be F. 

For the usual tail dependence case, Schmidt (2002) proved that when the 



radial random variable R has a regularly varying tail then X 1 and X 2 are 
tail dependent, and thus the tail order of the corresponding elliptical copula 
is k — 1. 

Example 1 (Student t copula) The radial random variable R for Student t 
distributions is a generalized inverse Gamma distribution such that R 2 follows 
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an inverse Gamma distribution with the shape and scale parameters being 
z//2, where v is the degree of freedom. It can be verified that Fr £ RV-j, (see 
Example 3 of Hua and Joe (2011a)). So the tail order for Student t copula 
1. 



IS K 



For univariate tail heaviness, one often use a concept referred to as Maxi- 
mum Domain of Attraction (MDA) of a univariate extreme value distribution. 
The following well known result characterizes the distributions that belong 
to the MDA of Gumbel, which is relevant to intermediate tail dependence 



of elliptical copulas. We refer to Embrechts et al. (1997) and de Haan and 



Ferreira (2006) for more details about MDA. For notation, A is the Gumbel 



extreme value distribution for maxima, and $ a is the Frechet extreme value 
distribution with parameter a > 0. 

Theorem 2 A random variable X with cdf F is said to belong to the Gumbel 
MDA (denoted as X e MDA(A) or F e MDA(A) ) if and only if there exists 
a positive auxiliary function a(-) such that 



.. F{x + ta{x)) _ t 
hm _ = e , 

x^oo F[x) 



t E 



(4) 



where a(-) can be chosen as a(x) = J°° F(t)/F(x)dt. 



For the case where R has lighter tails than any regularly varying tails, 
some asymptotic study has been conducted for elliptical distributions where 
R G MDA(A). We refer to |Hashorva| Q2007D, |Hashorva| (|2010[) and |Hashorva 



( 2008[ ) for relevant references. Now we are ready to present a result that is 
useful to find the tail order of a bivariate elliptical copula where the radial 
random variable R belongs to Gumbel MDA, and this result is cited from 



Hua (2012). Hashorva (2010) has a version of this result in Theorem 2.1 but 



the proof here is different. 

Proposition 3 Let C be the copula for an elliptical random vector X := 
(X 1 ,X 2 ) constructed as {J, and b e = ^2/(1 + g). If R e MDA(A), then 
the upper and lower tail orders of C is 



lim 



log(l-F fl (6 g r)) 
« hg{l -F R (r)) 



(5) 



provided that the limit exists. 
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Proof: Letting r := F 1 (1 — u) and b g = ^2/(1 + g), then by Example 6.2 



i) of Hashorva (2007), as u — > + and thus r — > oo, 



C(l -u, 1 -u) 

= P [X, >F-\l-u),X 2 >F-\1 

[1 - q 2 Y' 2 



U 



2?r(l - £) 2 



= P[Xx >r,X 2 >r] 
a(b e r)/r][l-F R (b e r)}, 



(6) 



where Fr is the cdf of R and a(-) is an auxiliary function of R with respect 
to the Gumbel MDA in the sense of @. As u —7- + , i.e., r — >■ oo, both 



a(b g r)/r — > (see Theorems 3.3.26 and 3.3.27 of Embrechts et al. (1997)) 



and 1 - F R {b Q r) -»■ 0. Let := 1/[1 - F R (V)], then G : K ^ K + 

is increasing and the conditi on of Q is equivalent to that G G T-varying 
with auxiliary function a(-) (de Haan, 1970, Definition 1.5.1). The inverse 



function of a T- varying function is a II- varying function ( |de Haan , 1974 

/ •> t / \ 1 1 c\ \' x r 1 1 (W ' I 1 Tn nrnfnrn / ^ 1 /*"" I I t rn rTrinrr Aeon in M* "i - +- on nn ~^ri linn 



Corollary 1.10). Therefore, G G IT- varying. Assuming that an auxiliary 



function of G 1 is a (-), by Lemma 1.2.9 of de Haan and Ferreira (2006), 



the auxiliary function Oq(-) of the n-varying function G~ l is slowly varying 
at oo. Moreover, a (t) = a(G _1 (t)) (de Haan, 1974, Corollary 1.10). So, 
a(x) = a (G(x)). Then in (|6]), 

a(b Q r)/r = a (G(b e r))/r = a (l/[l - F R (b e r)])/r, 

while 1 — F R (b e r) is rapidly varying in r at oo due to the fact that G is 
T-varying and any T- varying function is rapidly varying (de Haan, 1970 
Theorem 1.5.1). Therefore, 

1 - F R (b r) = 1-F R ^2/(l + g)F-\l - u )) 

dominates the tail behavior of (j6j) as u — Y 0, and thus determines the corre- 
sponding tail order of the elliptical copula. By the definition of tail order in 
Definition [T], we may also obtain the upper tail order by the following 



k = lim 

u->0+ 



log C(l - u, 1 - u) 
logw 



log(l-F fl (6 r)) 
r->oo log (1 — F(r)) 



lim 



By Example 6.2 (iii) of Hashorva (2007), as r — > oo, 

P[X a > r] = (1 + o(l))(27r)- 1 / 2 [a(r)/r] 1 / 2 [l - F R (r)]. 
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Due to the similar argument as before, 1 — Fr(t) dominates the tail behavior 
of P[Xi > r], as r — > oo. Therefore, we may write 



lim 



log(l-F fl (6 e r)) 
r^c log(l - F R (r)) 

which completes the proof. □ 
It is very convenient to apply this method to derive the tail order if we 
know the tail behavior of R, and R belongs to MDA of Gumbel. By Theorem 
3.1 of Hashorva (2007), this result can also be extended to multivariate cases. 
For o?-dimensional exchangeable elliptical copula, of which the off-diagonal 
entries of S are all g and the diagonals are all l's, the tail order is ^ where 
b e is replaced by b 8t d = \/d/[l + (d — l)g\. 



Example 2 (Bivariate symmetric Kotz type (Fang et al. 1990) copula) The 
density generator 



Kx N ~ 1 exp{-#r*}, (3, «£, iV > 0, 



where K is a normalizing constant. By Theorem 2.9 of Fang et al. (1990), 
the density function of R is = 2irxg(x 2 ) = 2Ktix 2N ~ 1 exp{— f3x 2 ^}. So, 

the survival function is 



POO 

F R {x)= 2Knt 2N - 1 exp{-(3t 2 ^}dt 

J X 



Kir 



p- N/ *w N/ *- 1 exp{-w}dw 



Kit 



■) incomplete Gamma function 



p- 1 !™- 2 * exp{-/3x 2§ }, x -> oo, 



where the asymptotic relation is referred to Section 6.5 of Abramowitz and" 
Stegun (1964). Then by (O, we can easily get that 



k = bf = [2/(1 + 

Therefore, the tail order for the symmetric Kotz type copula is k = [2/(1 + 
g)]€. Gaussian copula belongs to this class with £ = 1, so its tail order is 
2/(1 + g) which is consistent to Example 1 of Hua and Joe (2011b). 
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5 Archimedean copula 



An Archimedean copula C has the following typical form 

C(u U ...,Ud) = V'O'Vl); • • • 



(7) 



where if) is referred to as the generator of the Archimedean copula C, and 
if) needs to satisfy certain conditions (some papers or books also refer <p = 
ip~ l as the generator; e.g., Nelsen (2006)). Most of the commonly used 
Archimedean copulas correspond to if> being a Laplace Transform (LT) of a 
positive random variable X (Joe, 1997; Joe and Hu, 1996); that is, if)(s) := 



/ °° exp{— st}F x (dt), where F x is the cdf of X. In this case, if) is completely 
monotonic, if)(x) is decreasing in x, if)(0) = 1 and if) (oo) = 0. When if) 
is of such a LT, it can generate Archimedean copulas of any dimension. 
However, this condition is not necessary for a ci-dimensional Archimedean 
copula with d being given and finite. A necessary and sufficient condition for 



finite dimensional Archimedean copulas has been given in Malov (2001) and 



McNeil and Neslehova (2009) 



The above two sets of conditions on if) correspond to two types of stochas- 
tic representations of Archimedean copulas, respectively. One is the well- 



known frailty model (Marshall and Olkin, 1988 Oakes, 1989), and the other 



2009) 



is the recently studied scale mixture model (McNeil and Neslehova 
(or the resource sharing model in the sense of Genest et al. (2011)) for finite 
dimensional Archimedean copulas. In this chapter, we will discuss condi- 
tions that lead to intermediate tail dependence of Archimedean copulas for 
both stochastic representations. The former has been studied in the litera- 



ture (Hua and Joe, 2011b), and the latter is only studied for the lower tail 



(Larsson and Neslehova, 2011). In what follows, we will first review existing 



studies on intermediate tail dependence for Archimedean copulas, then we 
will present our findings of conditions that lead to upper intermediate tail 
dependence of Archimedean copulas through the scale mixture model. 



5.1 Resilience or frailty models 

Write Gj(u) := exp{— (0 < u < 1) for j = 1, . . . , d, then clearly G/s 
are identical univariate cdfs. Then, 



C(ui, ...,u d ) = 
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where Fh is the cdf of the resilience random variable H, and ip{ s ) = ^-ff( s ) = 
Jo°° e ~ s ^Fn {dQ- The mixture representation means that: there are random 
variables X\, . . . , Xd such that given H — £, they are conditionally indepen- 
dent with respective cdfs G{, . . . ,G d . Let Fj(l — u) := Gj(u), then F/s are 

also cdfs, and C(ui, . . . , Ud) = / °° YYj=i Fj(^~ u j) -^#(^0- Assume that Fj is 
the cdf of Yj for each j. Then we can also look at the copula C as the survival 
copula for the vector (Y"x, . . . , Yd) that are independent conditioning on the 
frailty random variable H; that is why we refer to this representation as a 



frailty model (Genest et al. 2011). Actually, from different perspectives, the 



same random variable H can be referred to as either a resilience or a frailty 



random variable. We refer to Marshall and Olkin (2007) for more discussion 



about the concepts of resilience and frailty. 

We now use the resilience model ^ to study how tail heaviness of H af- 
fects the tail behavior of the corresponding Archimedean copula, and we refer 
to Hua and Joe (2011b) for relevant studies. For any given (m, . 

,u d ) 



. ,u d ),0 < 
and thus a 



Ui < 1, a larger value of £ leads to a smaller value of C(ui, 
heavier right tail of H tends to generate stronger positive dependence in the 
upper tail. Sufficient conditions on the tail heaviness of H have been given 
in order to get an intermediate upper tail dependent Archimedean copula. 
If we use Mh '■= sup{m > : E[i/ m ] < oo} to describe the degree of tail 
heaviness of H, then under certain regularity conditions, k < Mu < k + 1 
with k G {1, . . . , d — 1} leads to upper intermediate tail dependence. The 



next result is presented in Hua and Joe (2011b). 



Proposition 4 Suppose ip is the LT of a positive random variable Y with 
k < M Y < k + 1 for some k e {1, . . . , d - I}, and \^ k \0) - ^ {k) (-)\ is 
regularly varying at + with the associated slowly varying function £ satisfying 
lim s _ s>0 + K s ) < 00 ■ Then the Archimedean copula has upper intermediate 
tail dependence, and the corresponding upper tail order is kjj = My. 

There are some upper intermediate tail dependent Archimedean copulas 
that have a simple form. One example is the Archimedean copula constructed 



by LT (Joe and Ma, 2000) 



exp{-w a }^/r(l + a- 1 ), 



< a < 1, 



and the upper tail order is Ku — 1 + a and the lower tail order is Kl = d a . 
Another one-parameter Archimedean copula that has a very flexible upper 
tail dependence structure is the following 
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Example 3 (Archimedean copula based on inverse Gamma LT (ACIG) Hua 
and Joe (2011b)) Let Y = X^ 1 and X follows Gamma(a, 1) with a > 0, and 
then My = ot and the LT of the inverse Gamma distributed Y is 



a) = —s a / 2 K a (2^~s), s > 0, a > 0, (9) 
r[a) 

where K a is the modified Bessel function of the second kind. For a G (0, +oo) 
noninteger, Ku{C^) = max{l, min{a, d}} and kl{C^) = \fd. 

The interesting property of ACIG is that, it captures a very wide range 
of upper tail dependence patterns by involving only one parameter. That 
is, when < a < 1, there is usual tail dependence in the upper tail; when 
1 < a < d, intermediate upper tail dependence is present; when d < a, the 
upper tail becomes tail quadrant independent (i.e., K — d). 

The condition of My in Proposition [4] being noninteger seems to be un- 
necessary. In the next subsection, the restriction will be relaxed when we 
study the tail behavior of Archimedean copulas through the scale mixture 
model. 



5.2 Scale mixture models 

An Archimedean copula can also be represented as the survival copula for a 
random vector (McNeil and Neslehova, 2009) 



X := {X x 



,X d )^Rx(S 



S d ) 



(10) 



where R and Si are independent for % = l,...,d, R is a positive random 
variable and (Si, . . . , S^) is uniformly distributed on the simplex {x e : 
Y2i x i = !}• The relationship between R in (10) and H in ([§]) is given in 
Proposition 1 of McNeil and Neslehova (2010); that is, 



R = Ed/ H, 

where E d is independent of H and is Erlang(rf) distributed (i.e., E d follows 
Gamma(<i, 1)). The Archimedean copula can be constructed as 



C^ d {u x ,...,u d ) :=ip{4) 1 {u 1 )-\ Vi) l {u d )), 



(11) 
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where the generator ip is the Williamson <i-transform of cdf Fr with -Fr(O) = 
(Williamson 1956); that is, 



(1 - s/r) d_1 F(dr), s€[0,oo) 



Note that Williamson ^-transform of a positive random variable can also lead 
to a generator ip that is completely monotonic (see Example [4]) . Throughout 



this subsection, assume (p(x) 
C$ td , then 



If the joint cdf of (Ui, . . . , U d ) is 



(X 1; ...,X d ) = (^(C/i), . . . , p(E7 d )) = x(S u ...,S d )] 



that is, for each margin Xj, P\Xj > x 
1 in 



Larsson and Neslehova 



where Y ~ Kumar aswamy( 



ip{x). From the proof of Theorem 

(2011), we know that X; = 

1 



BY for z = 1, . . . , d, 
[i-x) d -\x e [0,1]. 

We refer to Song and Gupta (1997) for more discussion of margins of a L p - 



, d—1); that is, F Y (x) 



norm uniform distribution. 



Recently, Larsson and Neslehova (2011) has studied the tail behavior 



of Archimedean copulas via the scale mixture representation. In Section 5 



of Larsson and Neslehova (2011), some tail dependence patterns have been 



derived except for the intermediate upper tail dependence case. We will fill 
the gap in this subsection. 

The lower intermediate tail dependence of bivariate Archimedean copulas 



has been studied in Proposition 7 of Larsson and Neslehova (2011), and the 



conditions needed on tp is essentially the same as Theorem 3.3 of Charpentier 



and Segers (2009). A more intuitive and fairly general pattern of ip has been 
considered in Hua and Joe ( |2011b[ ); that is, as s — > oo, ip(s) ~ T(s) = 
a\s q exp{— a>2S 1 ~P} and i[)'(s) ~ T'(s). In this case, if < /3 < 1, then 
has lower intermediate tail dependence with 1 < kl(C^) = d 1 " 13 < d. 



With survival copula for (10), the upper tail of R may influence the lower 



tail of the corresponding Archimedean copula. Using the notation of lower 
tail order k an d the relation (JTj), the re sult about the lower intermediate tail 
dependence in Larsson and Neslehova (2011) becomes the following 



Proposition 5 Let C be a d- dimensional Archimedean copula constructed 



as the survival function of the random vector in (10) and the associated ip in 



(11) is differentiable. Assume further that the radial part R in (10) satisfies 



R e MDA(A) with the auxiliary function a(-) of R satisfying a £ for 
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some < /3 < 1, and |ip holds, then Kl = d 1 13 and thus C has lower 
intermediate tail dependence. 



Based on (10), the lower tail of R, or equivalently, the upper tail of 1/R 
may affect the upper tail of the associated Archimedean copula. In what 
follows, we will prove that 1/R belonging to the MDA of Frechet (written as 
1/R G MDA($ Q )) may lead to upper intermediate tail dependence for the 
associated Archimedean copula; here the condition of the tail order being an 
integer is relaxed. 

Proposition 6 If 1/R G MDA($ a ) with k < a < k + 1, where k G 
{!,..., d — 1} is a positive integer, and ip is the Williamson d-transform 
of F R such that F R (0) = 0, and if k = d — 2 or d — 1, then further re- 
quire thatip is (d+l) -monotone and (d + 2) -monotone respectively, then the 
Archimedean copula constructed by ip as (11) has upper tail order kjj = a 
and thus upper intermediate tail dependence if 1 < a < d. 

Proof: First note that either ip being (d+ l)-monotone or (d + 2)-monotone 
can imply ip being d- monotone (McNeil and Neslehova, 2009[ ). So we can 
still apply the Willianmson <i-transform for the two cases where k = d — 2 or 
d — 1. By definition, the Williamson <i-transform of R is 



= / (1 - s/r) d - 1 F R (dr) = I >;(".*] (-s/ryF R (dr) 




(12) 

\ / / 

i=l 
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For any 1 < a < d, there exists a positive integer k G {1, 2, . . . , d — 1} 
such that 

k < a < k + 1. 

To study the upper tail for an Archimedean copula, we now investigate the 
behavior of the functions m,(s) in (12) for i = k, . . . , d — 1, as s — > + . 
Depending on the value of i, we consider the following cases: 
Case 1: a < k + 1 <i < d— 1. 

The condition l/R G MDA($ a ) implies that P[l/i2 > ■] e RV„ a , and 
i 7 /? = P[R < -] e RV a (0+). Write Fr(s) := s a £ fi (s), where £ B G RV (0+). 
Since y ^ F R {l/y) G RV Q , we have y i— >■ F R {l/y)y i - 1 G RVi_ Q _i. By 
Karamata's Theorem (e.g., Resnick (2007)), i > /c + 1 > a implies that 

.1/, 



F R {lly)tf- l dy 



1 



z — a 



and thus 



is 



l/s 



Fnil/yW^dy 



a 



Fr(s) 



a 



s~l R [s 



s -> + . 
(13) 



Case 2: 1 < % < a. The condition F R G RV Q (0 + ) implies that 



l/s 



F R (l/y)y l - 1 dy / EfiT*] < oo, as s ^ 0+ 
Therefore, mi{s) ~ E[i? _i ]s l as s — > + . 



Case 3: i = k = a. Let 4(s) := A; F R (l/y)y k 1 dy. Similar to the 



derivation of (13), by Karamata's theorem (e.g., Theorem 2.1 (a) of Resnick 
d2007| )), 4 G RV (0+), and hence, m fc G RV fc (0+). 

Then ^(s) = 1 + Y^Zi (V) (- 1 ) <m «( s ) implies that: 
if a is a positive non-integer, that is, k < a, then £j(s) — >■ E[i?~ l ] for 
z = 1, . . . , k, and 



i=l 
d-1 

+ E 



i=k+l 



d-1 

i 

d-1 

i 



i — a 



-s a £ R (s), s^0 + ; 



(14) 
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if a is a positive integer, that is, k = a, then 4( s ) — > E[R l ] for i 
1 k — 1, and 



(fc-l)Vl 



i + i{*>i}x yj 

i=i ^ ' 



-lYsH 



+ 



d- 1 

k 

d-l 



■l) l s*4(s) 



i=fc+l 



d- 1 



2 — OJ 



s ->• + . 



(15) 



Therefore, 1 < a implies that the map 1 — G RVi(0 + ). Write 1 — 
V>(s) = s£ (s), where £ G RV (0+). 

Assume that there exists a constant k and a slowly varying function 
4 G RV o (0 + ) such that C(l 1 - u) ~ u K £*(u) as u ^ + . If we can 

prove that k = a, get the expression of 4, and prove that such an £* is a 
slowly varying function, then the proof is finished. 

Let s := — u), then 1 — ip(s) = s£ (s) implies that 



1 = lim 



C(l -«,...,!-«) 



lim 

u->0+ 



i + E|=i(-i) i (J)V'(jV'- 1 (i-«))} 



[w K 4H] 
s)Ul-^s))}. 



The locally uniform convergence of 4 at + , together with 1 — ip(-) G RVi(0 + ) 
implies that, for any given t > 0, 



lim 



t(l-V(^)) 



t(^x(l-#)) 



lim 



lim 



1. 



>o+4(l-V(s)) «->o+ 4(1 --0(s)) «->o+ 4(«) 

Hence, 4(1 - V(0) e RV (0+). L et l(s) := lg(s)4(l - ip(s)), then £ G 
RVq(0 + ), due to Proposition 1.3.6 of Bingham et al. (1987). 



Applying the Monotone Density Theorem (Bingham et al. , 1987, Theorem 
1.7.2b) k — 1 times on (12) implies that ^^(0) < oo for i = 0, . . . , k — 1. 
Moreover, choosing u>j = 1 for each i in Lemma 2 of Hua and Joe (2011b) 
implies that Ej =1 (— l)"* (^) j l = for any positive integer i that is less than 
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d. Therefore, £* =1 (-l) j (J)jV W (0) = for i = 0, . . . , Jfe - 1. Define [x] := 
max{z integer; z < x}, and for any y > 0, y\ := y X (y — 1) x • • • x (y — [y]). 
By the l'Hopital's rule, 

1 = 1 . zU-^O^ 1 ^^) (m] 

s ™+ [k\/(k — k + l)\]s K ~ k+1 £(s) ' 1 ' 

If a is a non-integer, that is, if k < a < k+1, then applying the Monotone 
Density Theorem two more times for (16) with respect to (14) leads to 

d-l 



(-l)*+V (fc+1) (s) 



i=k+l 



d- 1 

i 



2 — a 



[a\/(a-k-l)\]s a 



-fc-i. 



(17) 



lim 



lim 



[k!/(k- Jfe)!]s*-*€(s) 



E 



•IV 



-k-l Id 



[«!/(«-*-!)!]* 



.re— fc — 1 / 



(18) 



Combining (17) and (18) leads to k = a, and then 

d 



j'=i 



d-l 

E 

i=jfe+l 



d-l 

i 



-V, 



a 



which is a slowly varying function. Thus, the slowly varying function £* can 
be chosen accordingly, which proves the case where a is a positive non-integer. 

If a is a positive integer, that is, a — k, then applying the Monotone 
Density Theorem one more time for ( 16 ) with respect to ( 15 ) leads to k = k = 



a, and similarly the slowly varying function 
□ 



can be obtained accordingly. 



Example 4 (ACIG copula) For the ACIG copula studied in Example 4 
of Hua and Joe ( 2011b[ ), H^ 1 ~ Gamma(tt, 1) for a > 0. Therefore, R has 
the same distribution as the product of two independent Gamma random 
variables with scale parameter 1 and respective shape parameters d and a. 
The product follows a K-distribution and the density function is 

f R (x;d,a) = —J—x^ 2 - 1 ^^), 



T(d)T(a) 

x G [0, oo); d positive integer; a > 0, 



(19) 
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where K is the modified Bessel function of the second kind. We refer to 



Jakeman and Pusey (11978) for the reference of density functions of a Re- 



distribution. Because the left tail behavior of R affects the upper tail depen- 



dence pattern of the corresponding Archimedean copula derived from (10), 



we need to study the behavior of (19) at 0. If a is not an integer, then 



K d _ a (s) ~ * (r(d - a)(s/2r- d + T(a - d)(s/2) d - a ) , s -+ 0^ 



Therefore, 



. F(d — a) n ■> Via — d) ^ -, 
f R (x;d,a) ~ T g + ^7-^7 V " 1 , a; -XT 



T(d)r(a) 



T(d)T(a) 



If 1 < a < d, then the term x a 1 dominates the tail behavior at 0; this is 
upper intermediate tail dependence case. When d — a is a positive integer, 



then by Abramowitz and Stegun (1964), 



d-a-l 



K d _ a (s)~l(s/2) a - d 



k=0 



■l) k {d-a-k-l)\ 
jfc! 



8/2) 



2k 



{d - a - 1)l (s/2T- d , s^O- 



Therefore, 



(d — ot — l)\„-i T(d — a) n i 
f R (x; d, a) ~ V / x"- 1 = { x a ~\ x^Q- 



T(d)T(a) 



T(d)T(a) 



which is the same as a being a non-integer. Combining these two cases, 
1 < a < d implies upper intermediate tail dependence of the ACIG copula. 



Proposition [6] is actually very useful to guide us to find a scale mix- 
ture random vector whose survival copula is an upper intermediate tail 
dependent Archimedean copula. For example, let R follow the positive 
Weibull distribution; that is, Fr(x) = 1 — exp{— x a } ~ x a ,x — > + with 
1 < a < d. Then Proposition [6] implies that kjj = a and thus the correspond- 
ing Archimedean copula has upper intermediate tail dependence. Another 
example of Archimedean copula that has upper intermediate tail dependence 
is presented in the following. 
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Example 5 (Dagum-simplex mixture) The Dagum distribution is also 
referred to as an inverse Burr distribution, and it is a special case of gener- 
alized beta distribution of the second kind (e.g., Kleiber and Kotz ( 2003[ )). 



Let the cdf of the radial random variable R be Dagum, then 
F R (x) = [1 + {x/a)- a ]~\ x>0, a,p,<r>0. 

We choose a = 1 for simulation as the scale parameter does not affect the 
associated Archimedean copula. It can be derived that F R 6 RV Q , / 3(0 + ). By 
Proposition [6j if 1 < a/3 < 2, then the copula C(u,v) := ■0(V>~ 1 ('u) +ip~ 1 (v)) 
should have upper tail order kjj = oif3. The simulated scatter plots are 
illustrated in Figure [TJ where the left plot is for uniform margins, and the 
right plot is for standard normal margins. The sample size was 2000 for the 
simulations. 



Figure 1: Simulation of Dagum-simplex copula 




In Figure [TJ the upper tail order is Ky = 1.08, which belongs to upper 
intermediate tail dependence. For the lower tail, Fr(x) = 1 — (1 + x~ a )~ 13 , 



and F R (x) 



-a/3(l + x- a ) 



-J8-1 



X~ 



—a(3x 



as x — > oo. So F R e 
RV_ Q . By Corollary 2 of Larsson and Neslehova (2011), there is lower tail 
dependence and the tail order parameter Xl — 2~ Q . 



20 



6 Remark and future work 



The notion of tail order provides a quantity to evaluate the degree of de- 
pendence in the tails of joint distributions, especially when intermediate tail 
dependence appears. We first review fundamental concepts and existing re- 
sults of intermediate tail dependence. Throughout the review, some new 
properties of intermediate tail dependence have been given to supplement 
existing results. The new results mainly consist of an easy way to derive 
the tail order of a bivariate intermediate tail dependent elliptical copula, and 
the study of intermediate upper tail dependence for Archimedean copulas 
constructed from a scale mixture model. 

Proposition [6] is helpful for constructing an upper intermediate tail depen- 
dent Archimedean copula. However, the scale mixture approach can often 
only give us a simple way to simulate desired tail dependence structures, but 
not a simple closed-form parametric copula family. So, how to apply the scale 
mixture model to provide various desirable models and to make statistical 
inference efficiently will be a very interesting topic for future research. 

As mentioned earlier in Section [TJ there is a link between tail orders 
and multivariate stochastic orders. Future research also includes whether 



certain forms of multivariate stochastic orders (Shaked and Shanthikumar 



2007, Chapter 6) can be adapted to tail forms for comparing the strength of 



dependence in the tails of copulas. 
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